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Surface results are compared to those of a Born-Mayer/Morse 


SBetencszal function. Calculations show that the effects of 
vacancies exceed those of adatoms, requiring a surface 
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Change in the yield. The choice of the potential function 
affects the sputtering properties to the same degree as the 
surface defects, and tends to cast some doubt on 


quantitative results from this type of simulation. 
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a INTRODUCTION 


Ars Whewl DPS SPUPPERING? 

The ejection of surface atoms from solid surfaces under 
bombardment by energetic particles is known as sputtering. 
The particles may be i10neg, neutral atoms, neutrons, 
electrons or photons. Ejection aes causedemby collisions 
between the incoming particles and atoms in the selvage, 
id.e., the surtace layers of the solid. The processes 
involved in sputtering are, in principle, similar to those 
causing radiation damage in the bulk of a solid. These 
processes occur far trom thermal equilibrium, which means 
that sputtering is different from evaporation, which does 
occur at thermal equilibrium. 

Generally, an incoming particle will collide with the 
atoms of the solid, thereby transferring energy to the 
atomac nuclei. If the energy transferred to an atomic 
nucleus is greater than the binding energy at the lattice 
site, a primary knock-on atom (PKA) is created. The PKA 
Will collide with other Seo ee atoms distributing the 
energy by way of a collision cascade. An atom is sputtered 


if the energy transferred to it has a component normal to 


iat 


the surface which is larger than the surface potential 
energy barrier’. 

One measure of sputtering is the sputtering yield, Y, 
defined as the mean number of atoms removed per incident 
particle. In counting the removed atoms, only those of the 
solid are included, while reflected 16nsPon anci@emt avens 
that have been reemitted are not taken into account. The 
particles removed by eputtering are emitted with a broad 
distribution of energy, in different excitation and charge 
states, and at all exit angles. Differential yields are 
defined to describe the variation of the yield with these 


parameters. 


Be APPLICATIONS OF SPUTTERING 
Sputtering was long regarded as an undesirable effect 
which: 1) destroys cathodes and grids in gas discharge 
tubes or ion sources, a2) contaminates a plasma and 
Surrounding walls, making it a major research area in 
fusion technology, anes saba: 12) causes the destruction of 
diaphragms and targets in accelerators and in high-voltage 
electron microscopes. However, sputtering is used today 
for many applications and has become an indispensible 
process in modern technology. 
, The potential energy of a surface layer atom is 
normally referred to as the surface binding energy (SBE). 
The SBE, however, only has meaning in an undamaged lattice 


and does not define the energy required to remove a surface 
atcen once the collision process has begun. 


ee 


spettering allowe a controlled removal of tightly bound 
Burface layers on a nearly atomic scale and possible 


submicron spacial resolution if a well focused or rastered 


beam is used. The atoms removed can be analyzed in a mass 
spectrometer, which givee Toor neaeltOn asour suritace 
concentration, or can provide a depth profile if the 


Ssurtace is continuously bombarded. 

Both the removal ot atoms from a surface and the flux 
of atome leaving the surface are successfully applied, for 
example: in eputter ion sources, to obtain atomically clean 
eurtaces; in micromachining; ana Serteoremaesthn ~ orotiling of 
thin films and surfaces. One ot the largest applications 
of eee is the deposition of thin films ona large 
variety of substrates. These filme may be several square 
meters in area, or may be extremely small as in 


microelectronics. 


e HIST@RICAL OVERVIEW 
Sputtering has been investigated in most detail for 


bombardment of monatomic solids with mercury ions, with 


noble gas ions and with hydrogen ions. Sputtering yields 
may lie between O and 10* atoms per incident ion, but are 
Motcally one to five. The yields depend on: 1) the 


incident particle’s energy, mass and angle ot incidence; 


ii) the mass of the target atoms; iii) the range of crystal 
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order® and the crystal orientation of the solid and iv) the 
surface einaiee energies of the target. The yield is, 
however, nearly independent of the temperature [Refs. 1, 2, 
od Below a threshold energy, which is about 20-40 eV for 
normally incident ions, no sputtering takes place. Above 
this threshold the yield increases with incident ion energy 
and reachs a broad maximum in the energy region of 3-50 
keV. At higher energies the sputtering yield decreases. 
This decrease is related to the larger penetration of the 
ions into the solid and the lower energy deposition into 
the surface layers. 

Several versions of the theory of sputtering currently 
exist in the literature, each ot which exhibits the 
propensities of their authors. Each version has some 
experimental support, but certain approaches are in better 


repute than others. This review makes no attempt to 


Uf 


rationalize the various approaches, but. raetherea@it 4 
concerned with providing a historical overview of the 
development of the theories. 

Sputtering was first recorded in teS2) by Greve ehewe: 
4, 3] who noticed the disintegration of cathodes in glow- 
discharge tubes with subsequent deposition of the material 


On the surrounding glass walls. This became known as 


. The range of crystal order refers to the distance 
over which the crystal maintains the regular, repetetive 
Qridlike arrangement of atoms, known as a lattice. 
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‘aath@edem aputtemdng" in thew literature. It took almost 
tifty years before Goldstein (CRef.6] presented compelling 
evidence that the sputtering effect was caused by positive 
ions of the discharge striking the cathode. 

Untortunately, early investigations of sputtering were 
wasted before it was realized that the pressure surrounding 
the target was a critical parameter and had to be low 
enough to allow the sputtered atoms to escape the surface. 
Many anomalies in early experimental results can be 
accounted for by the simple fact that the mean free path of 
the sputtered particles was less than the distance to some 
“esltlector surface", which allowed the liberated particles 
to diffuse back to the surface after collision with a gaé 
atom. Penning and Moubis (Ref. 71] were the first to 
conclusively demonstrate the effect of ambient pressure, 
even though the view that pressure was an important factor 
had been held for a few years prior to their publication. 

The work of Guntherschultze and Meyer (CRefs. 8, 9] is 
one of the few early experiments preceding the findings of 
Penning and Moubis which satisfied the conditione for a 
reproducible sputtering determination.- Guntherschultze and 
Meyer operated their discharge tube at relatively high 
vacuum (for that era) with a copper target and took the 
precautions of removing initial layers of the target 


material. 


naka 


The elimination of adsorbed gases and surface oxide 
layers as a critical factor in deterring ne weer u- a 
wae further described by Ariiov, et “alo, tn soem net. Se 


Yonts and Harrison also presented evidence (Ref. i1] that 


surface recontamination from  beackoqroaund gases was a 
Signitiecant —reeccor aka) quantitative sputtering yield 
measurements. Criteria were theretore developed tor 


conducting "clean" experiments. For a single crystal metal 
this means the removal of surface and near-surface damage 
and impurities by mechanical abrasion and polishing, with 
thermal annealing and "sputter cleaning" used to remove any 
residual damage. It is possible to go to greater lengths 
to ensure damage- and contaminant-free surfaces, oe 
example, using low energy heavy ion bombardment to remove 
surface oxide (conducted in ultrahigh vacuum) followed by 
thermal annealing to remove surtace damage caused by the 
ion bombardment. Intormationvaboute the atonicrst hue euge see 
the surface can then be obtained by low energy electron 
diffraction, Seep field ion microscopy of the new 
tunneling electron microscopy. Subsequent sputtering 
experiments could then be conducted in an ultrahigh vacuum 
environment (less than 10-'® Torr) to preclude adsorption 
of gas molecules on the surface. 

The concept of an individual sputtering event on an 


atomic ecale wage firet recognized and analyzed by Stark. 
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He originated the so-called hot-spot model? of sputtering 
CRet. 12], and subsequently, a collision theory viewing 
sputtering as a sequence of binary collision events 
feeeteaced by One bembarding dion at a time (Ref. 13]. 
Metatzing the conservation laws of the theory of elastic 
Semis lons and the Goncept of collision cross sections in 
his analysis, Stark correctly interpreted the observed 
energy dependence of the sputtering yield of hydrogen ions 
bombarding a metallic target. 

Subsequent investigators took a ditferent stand on the 
issue. Atter Kingdon and Langmuir successfully applied 
Seerk’ Ss eollieion theory PRet. 14], the erroneous 
impression arose that the collision theory would imply that 
Sputtering was a single-collision process resulting ina 
Sememagly peaked angular distribution of the ft hUeeot 
sputtered particles. Only after the demonstration by 
mem@er tChefe. 15, 16] of crystal structure ettects in the 
Ga Gt sputtered particles did it become evident that 
local evaporation alone could hardly explain the sputtering 
phenomena. Wehner discovered that ae angular 
distributions of the atoms ejected from single crystal 


targete showed maximum intensities in eS Sc Ease = 


a The hot-spot model was considered to be the 
evaporation of target material from a microscopically small 
region with a high local temperature due to individual ion 
impact. 


1/7 


corresponding to the more closely spaced atomic rows 
passing through the target surfaces. 

The discovery of such ejection patterns or "Wehner 
Spote" revived the interest in the collision theory of 
Sputtering. Independently, Keywell (Refs. 17, 18] made a 
first attempt to formulate Stark’s multiple-collision model 
in terms familiar from neutron transport theory. Keywell’s 
work, as well as subsequent calculations by Harrison (Ret. 
19] were important steps in the sense that probability 
concepts as expressed by collision cross sections finally 
made their entrance into sputtering theory. 

The theory ot collision cascades has been further 
advanced by the work of Leibfried [Ret. 20], Lindhard, et 
al. CRef. 211], Dederichs (Ref. Bey, Robinson (Ref. 231], 
Thompson (Ref. 24], and Sigmund (Reis; 23f26), wheenac 
alsa collected the available knowledge in a trangport 
theory of ‘sputtering bret. 277 In an amorphous solid the 
collision cascades have been described by a Boltzmann 
transport equation. Sigmund was able to obtain first order 
asymptotic solutions for the eee: cascade regime’. 
However, only polycrystalline solids with randomly oriented 


crystallites can be approximated by an amorphous solid. 


; The linear cascade regime applies to bombardment 
With ions ot medium to large atomic number in the keV 
energy region. Here large collision cascades can develop, 
however, moving target atoms are consideréd to collide only 
With target atoms at rest. 


18 


Despite its shortcomings, this theory hag become successful 
as a reference standard for sputtering yield measurements. 
It has also been a reasonable starting point for further 
refinements of the theory. During the rapid development of 
the theory of collision cascades, the hot-spot (or "spike") 
theory has received little attention despite repeated 
claims of a thermal component observed in energy spectra 
Ot sputtered atoms [CRet. 28]. 

Macnee periledg trom 1955 to i965, the dominant effort in 


the study of sputtering on metals was actually spent in the 
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Pees tigatiton of crystal lattice effects. These ettfect 
fe inte twO distinct groups, the channeling theory of the 
Sputtering yield and the focusing collision theory of 
ejection patterns. The channeling theory accounte for the 


aeeulear variations of the yield cloge to the principle 


crystal axes and planes, but is not very successful for 
more general directions. The tecusing =ecollision eheary 


Gives a qualitative description of the ejection pattern 
spots from cubic metals at high bombarding energies, but 
cannot successfully interpret low energy irradiation 
resutts or observations on hexagonal or other low symmetry 
materials. The discovery of egpot patterns by Wehner (Refs. 
15, 16) inepired Silsbee’s identification (CRef. 29) of the 
focusing collision sequence as a means’ for long-range 
transport of momentum in crystals at low energy. An 


alternative to the focusing theory of the ejection pattern 


ie, 


epote wae propozed by Lehmann and Sigmund (Retr. 30). ire 
model stipulates that the target surtace have an ordered 
etructure, but not necessarily long straight rows of atoms 
intersecting the surface. This model emphasizes the role 
of surface structure and of the surface binding energy in 


monocrystal sputtering. 


De COMPUTER SNBDEEING 

In the past, sputtering theory needed only to predict 
the sputtering yield as a function of energy for a given 
incident particle and polycrystalline target of a single 
metal. The experimental environment led the theorists to 
Statistical theories involving ensemble averages over 
distribution ftuncewons, ultimately culminating in the work 
of Thompson (Ref. 24) and Sigmund (Ref. 25] which predict 
the yield reasonably well. However, the chief limitation 
ot these theories is their inability to deal with 
experimental results attributable to the regular structure 
in a Single crystal target. As is the case with most 
theory, the statistical sputtering theories are also 
constrained by the lack of analytic tractability. 

With the advent of high speed computing capabilities it 
became practical to include the full three-dimensional 
structures of crystals in the theoretical models of 
sputtering. The investigations of Vineyard and his 
colleagues (Ref. 313] demonstrated the power of computer 
Simulations in isolating and understanding the elementary 
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Peeeemsee OL detcet production and migration in radiation 
damage. Such techniques were first applied to monocrystal 
sputtering problems by Harrison and his co-workers (Ref. 
Beil. They proposed ejection mechanisms to account for the 
Wehner spots which differ from both the focusing collision 
theory and the model of Lehmann ane —oLomunG. 


These computer eimulations can be divided into two 


categories. “Historically, the binary collision (BC) 
Simulation models were developed first. The binary 


collision simulation is known as an "event store" model. 
The program tollows 4 cascade of moving particles 


sequentially, moving only one particle at a time, but 


remembering events, es > eolliscionse, as it progresses. 
The modele usually invoke the binary Ged & Lom 
approximation (Ref. Soe This assumes that each particle 


interacts with only one other particle at any time, with 
the other particle usually considered to be stationary. 
Thus these models are inherently linear calculations and 
theretore they have many of the same restrictions as the 
linear i ee calculations. 

Initially, these models simply mimicked the statistical 
theories, then Robinson and VDén included the target crystal 
etructure (CRef. 331] and ettfectively discovered channeling. 
The channeling concept fed back into sputtering theory and 
Onterderlinden CReft. 34] introduced the "channeling" 


sputtering model, which satisfactorily explained the 


veal 


relative magnitude of the egputtering yield from single 
crystal target surfaces at low keV ion energies. Since 
their inception, the BC programs, MARLOWE (Ref. 351], een 
[Ref. 36 J and MORLAY C[Ref. EVGA Me have assumed that the 
physical model was well understood and that the simulation 


Should be designed to compute results from it efficiently. 


it 
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Each program has its particular strengths and eftficienci 


because each was designed to address specitic criteria. 


The second method of simulation is the multiple 
interaction Cri Simulation and is based on molecular 
dynamics. MI simulation was motivated by the work of 


Gibson, Goland, Milgram and Vineyard (Ref. 31], with tie 
proponents et this model having concentrated om 
understanding the fundamentals of ejection processes. 

This type of simulation is known as an "timestep" mode 


and involves solving Newton’s laws of motion, often 


expressed in Hamiltonian form, for many particles. These 


Calculations proceed sequentially in time, following the 
collisions of many particles simultaneously. An advantage 
of the model is that there are no built in constraints as 
to the behavior of a collision cascade and, given enough 
cpu time and memory, will accurately te lteweer fone ndveod 
cascade to its ultimate thermalization. 

A new development is that of a hybrid between the two 
typee ot codes by Harrison, et. al. (Ref. 38]. The goal of 


which is to have a Simulation which retains the "exact" 


treatment of the time step model in the surface region and 


minimizes computation time for the remainder ef the 
cascade. 

As the Simulations have become increasingly more 
sophisticated, attention has Burned to Squad yameg 
Semiconductor materials. These materials present a 
Gmeeincc praeblem due to their nonspherical interatomic 
potentials. Computer modelling of silicon and other 


semiconductor materials may soon be a reality due to the 
pioneering work of B. J. Garrison in the use of embedded 


Beeemtcial tunctionse*. 


= Private communications between Professor 068. ay, 
Garrison of Penngylvania State University and Professor 
eee eadiaricem, Jr. af the Naval Postgraduate School. 
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i ines OBIE GRIVES 


Numerous experimental and theoretical investigations of 
the energy and angular distributions of alems and emeleeiwiee 
ejected from single crystal surfaces due to ion bombardment 
have been made over the years. Wehner was the first to 
observe spots, Or anisotropic angular distribuevenow eee. 
LS) fee His method of detection, however, did not allow 
him to distinguish the mass or energy of the ejected 
Species. In more recent times energy and angular 
distributions of ions have been measured by angle-resolved 
Secondary Ion Mass Spectroscopy (SIMS) CRefits. 39, 40> 4s 
42]. The theoretical investigations, mainly computer 
Simulations, have been used to predict the distribution of 
neutral species C[Refs. 43, 44, 495, 46]. These molecular 


dynamics calculations of the ion impact event have helped 


Loni) formulate the mechanisms ot particle emission from 
Single crystal surfaces, ii) provide an atomistic picture 


of the damage created within the material itself, and iii) 
unravel the origin of angular anisotropies. 

Experimental measurements of particle trajectories have 
not been possible at the level necessary to make direct 
comparisons to theoretical predicuitens: However, 
qualitatively and semi-quantitatively favorable comparisons 


have been made between the predicted nevtral distri put rame 
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and the measured ion distributions C[Refs. 40, 41, 42]. 


Rigorous quantitative comparisons are difficult since the 


Mica tton probability is not known. The ionization 
probability is required since the ions experience a 


ditferent interaction potential than the neutrals as they 
leave the surface. In particular, it would be advantageous 


to be able to directly compare experimental and theoretical 


0 
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distributions in order to use the measured distributions to 
Sees eruectural inteormatioh about the surface and to 
improve the theoretical model s0 as to gain a better 
understanding or the scattering process. 

With the recent advent of the energy- and angle- 
resolved neutral-particle Clams tector by N. “vyinogred, 
jen a | « CRef. 471], distributions for rhodium atoms ejected 
from clean and oxygen covered Rh(111} surfaces due to ion 
bombardment have been obtained [Refs. Hepes SO Is Energy 
and angular aistriputions reported by this method 
Semestrerute the first reported neutral-particle ejection 
distributions for material desorbing from genuinely clean 
well-defined surfaces. Surface cleanliness and structure 
were monitored by low energy electron diffraction (LEED) 
and by Auger electron spectroscopy (AES) as well as by the 
Reoredueri mai lity of the EARN measurements. 

iaee results of computer simulations of argon ion 
bombardment ot Slot st ab seal 3; Sti do net correlate 


Buameitatively with the reported EARN results. Thais. baek 


Ze 


ort Bee has provided the impetus to investigate 
possible mechanisms that could be the cause of the 
divergence between the simulated and experimental results. 
Qne of the assumptions that is commonly made in computer 
simulations is that the target surface is completely free 
of defects, i.¢., atomically emoocth. The degree oft surface 
detail -previded hy | EEED Vande ear ae EARN experiments 
does not preclude the presence of small scale topographical 
defects in the form of adatoms or vacancies, The only 
method that is currently available that could provide more 
precise characterization of the surface detail Wast hole 
disrupting the surface itself is scanning tunmaeling 
microsceapy (SIM) [Rets. sa) eed: 

This study attempts to expand the scope of variables 
that are normally considered in computer simulations by 
investigating the Si8ee OF small scale surface defects, in 
the form of single and multiple adatome and vacancies, on 
quantities such as the yield and the angle and energy 
distributions of SspuEEteread partie le—. The “Rhiggiel 1 vee ee 
System was chosen for investigation because the EARN 
experiments have set a new standard as being the most 
accurate, reproducible results obtained for sputtering to 
date, and provide the best baseline of experimental 
informa. aon. No attempt is made at this time to directly 


compare the results of the classical dynamics model to 
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experimental results, this being left for future studies 


when more E@xperimental data becomes available. 


di Bee h THE N@bEE 


A. MULTIPLE INTERACTION “Hio Strib 

The MI simulation used to model the ion bombardment 
events is QDYN, developed by Professor Don E. Harrisen, er. 
at the Naval Postgraduate School. A detailed description 
of the computer simulation method, the mathematical model 
it is based on, and the integration scheme have been 
published previously by Harrison, Site. aia (Refs. (aay oo- 
54]. 

Briefly, QDYN is a sputtering simulation which proceeds 
through a series of timesteps, each of which consists of 
the following: i) summation of the pairwise forces for each 
atom; ii) calculation of new velocities and positions at 
the end of the timestep; iii) movement of atoms to their 
new positions; and iv) test of energy conservation. 

QDYN computes the interaction forces by performing a 
numerical solution of Hamilton’s classical equations of 
MOLLon, using a predictor-corrector integration scheme 
which bases the time increment of each step on the fastest 
MOVIingG partic le, This differs from the tform oi mest other 
numerical integrators because position and velocity are 
calculated for the same instant Sse time. More 
eophisticated integration schemes do not improve the 


Calculations and only increase the computation costs. 
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The system consists of an incoming argon ion and a 
Single-crystal rhodium target containing up to 1854 atoms. 
meen particie is characterized by i) its mass, ii) the 
ence laws by which it interacts with other particles, iii) 
mes Position, and iv) its velocity. Hace tne wncoming Lon 
strikes the target at a specified impact point the 


positions and momenta ot the ion and all the atoms ain th 


iT 


Single-crystal target are allowed to develop in time. This 
Savewmlatioqn, called a trajectory, is a simulation of a 
Single collision cascade. The collision cascade is 
terminated when the momentum has dissipated through the 
lattice and no more atoms can be ejected. Pa taiesumetaltiec 
rhodium system, tiie corresponds te wea Maximum particle 
energy, in the solid, Or 220 ey. Finite lattice 
temperature effects and energy losses due to electronic 
excitations have not been included in this calculation. 

Two ancillary programs were used to investigate each 
Simulation. ANMOL analyzes the final positions and momenta 
to determine yields, the energy distribution, possible 
multimer formation, and the angular distribution of the 
ejected atoms, all ot which are experimentally measurable 
quantities. In addition, information about atoms ejected 
per single ton (AST) and ejection times is provided. 
Sranh ica 1 representations et the distributions and 
Quantities mentioned above were provided by ANPLOT. An 


alternate output available from QDYN allows sufficient 


Zo 


information to be printed out so that the col bstem, cagcade 


can be traced. 


Be PHYSICAL CONSIDERA UONS 
ik. Lattice Orientation and Constants 

A Rh{iil} orientation was chosen to match the 
conditions of the energy- and angle-resolved neéutral- 
particle (EARN) experiments (Refs. 48, 49, 50]. Rhodium is 
a face centered cubic (fcc) metal with a lattice constant 
(Ga) Gi 3, S042. Other physical constants of rhodium used in 
the development af a Rh-Rh potential function, incluce: Ene 
bulk Sheat of [eatonitzs ean (35,76 eV) and the bulk 
compressibility ©3.625%1057#en-7 <2). 

The basic distance unit defined for this simulation 
is the lattice unit (LU), which is one half the latttee 


constant for cubic materials. 


Pa Lattice Size and Containment Criteria 

To model the S keV bombardment of the EARN 
experiments and describe most of the ejection events would 
require a target size that is computationally not feasible. 
In order to gain the most complete information about the 
Rh{ilid/Ar+ system this study conducted its computations at 
ion bombardment energies of S00 eV and 2 keV. These 
energies and their corresponding lattice size were chosen 


in -e@rder “vane optimize containment of the collision 


20) 


cascade within the lattice, ie) maximize the number of 
atoms sputtered to ensure statistical confidence in the 
results, and iii) minimize the computational time necessary 
to complete a simulation. 

The inability of the Jattice of target atoms to 
contain the entire cascade created by the incident ion is 
Salled a tailure at "containment". her mines ple; there 1s 


no problem for a moderate-siz@e® camputer to store the 


information for even tens of thousands of atoms, ae. aL ye 
Sesame elee, running time considerations Gea, Meaney) bin © 
the target to roughly 2000 atoms. Cascades created by keV 


ions are known to involve many more than 2000 atoms, so the 


problem seems insoluble. Atter long experience with these 
computations, Harrison has defined a more practical 
@etaniti@on of containment [Ref. Sak: "If an increase in 


target size does not change the results of the computation, 
the trajectory is etfectively contained for the purposes of 
Picta Gomputation"™. 

ines praecice, for any observable, one finds that 
containment is approached asymptotically as the target size 
increases. Harrison also states that: 


As one might anticipate from experimental experience 
with absolute values, yield provides the most severe 


test af containment; 20 absolute yield containment 
guarantees eéverything else. As the target size 
increases, all othér global results reach constant 


values for much smaller targets than those required ta 
produce stable yield values. 
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This fact is Bern out by Figure 1. wite wea 
graphical summary of numerous calculations done on clean 
Rh{iii>} at various ion energies and lattice sizes. These 
calculations provide the baseline information on clean 
Rh{1ii?t as well as acting as verification of containment. 

The Po laws me lattice sizes (see Table 1) 
correspond to the information provided in Figure l. The 
lattice dimensions are given in terms of the number of 
lattice planes in the x, v; and 29 directroens. Figure 2 
provides a graphical description of the lattice plane 
numbering system. 

TABLE 1 
LATTICE DINENS LANS 


Lattice Dimension # Rh Atoms in the Lattice 


Lea Sec S341 
Lea Bee 648 
2k eS ees 683 
21. % Goes eco 
Zon Ge ee Dea6 
2a < -"Gaex  ee 276 
27 “SO eee 1540 
Zo x Oe 17 OL 


The lattice sizes chosen for the computations were: 


Lon Sener ay Vatlice S122 
900 eV 1S: Seen 
2000 eV Zo Xmen, 


Even though appreciable movement of the edge atoms 
was observed prior to the last atom leaving the surface, 


the signiticant features of the dynamics were not altered 


oe 


when a larger lattice size was used. All simulations were 
run on an IBM 3033 at the Naval Postgraduate School with 
computation times on the order of 1 CPU hour for the 500 eV 
Simulations, and 3 CPU hours for the 2 keV simulations. 

For the 18x5x12 lattice there are 108 atoms in the 
aaeoe ayer and for the 25x6x17 lattice there are 213 
atoms. To aid an describing the simulation events all 
atoms (and ions) in the simulation are numbered, with the 
10on being number 1 and the lattice atoms numbered beginning 
with number 2 at the origin and progressing upward across 
the rows, one layer at a time. Figures 3 and 4 delineate 
eae membering system for the tfirst layer of each lattices 
system. 

Another containment consideration deals with those 
atoms that leave the microcrystallite through the bottom 
and sides. The Simulation program considers these atoms to 
pemelose to the collision process, but continues to track 
them for energy conservation purposes until the model is 
EoGmoma ted. This is justified for those atoms that depart 
Premlattice through the bottom, Since almost all of the 
e@gerisiommevents that contribute to sputtering occur in the 
fi~et foumelayers of the lattice. The atoms that are 
ejected from the sides of the crystal are ignored in order 
to improve computational efficiency. Hater asan, through 


long experience with this model has determined that the 


eye: 


deletion of these atoms has little impact on the number of 
atoms sputtered, altering the results by about O.1%. 

Alternate approaches to those atoms lost through 
the sides of the microcrystal would be: i) to reflect them 
at the side boundaries, or iil) to reintroduce them on the 
apposite side of the lattiege from whieh they exited; with 
the same momenta. 

ae Statistical Ensemble 

An ensemble average is an average of any property 
over a group of similar systems. Simulations produce 
ensemble] averages, where each trajectory is a sample from 
an ensemble of all possible trajectories realizable trom a 
particular atom ejection. For a monatomic ion, each impact 
point on the target surface leads to a different 
trajectory. Using point and line symmetry transformations, 
one establishes a representative area in the target, which 
is an elementary symmetry zone ot the surface of that 
lattice orientation. An ensemble of trajectories is 
developed by uniform sampling otf the representative area. 

For our purposes a rectangular representative area, 
of dimensions Q.707 LU By 1.225 LU, was chosen, with its 
origin placed at the center of a target atom (Figs. 3 and 
4). Uniform sampling of the representative area was 
accomplished by desiqnating 300 impact pointS within the 


rectangle (Pig wa, 
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In order to reduce computation time, the necessity 
Seeing a Unitorm distribution of 300 trajectories was 
investigated. The impact points were divided into two 
groups of 150 points each (indicated by x’s and +’s in 
rose. 2) and the results of separate simulations of 500 eV 


Ar* on Clean Rh(1il)} were compared (Table 2). 


TABEE: 2 
S60 "eV ENSEMBLE RESUS 


Number of trajectories 300 PoaoeCs so) a0 (6+ s 
*# Rh atems sputtered =| pe. 457 445 
ASI Sa 3,708 eee 7 
Average energy of 

eputtered Rh (eV) Soue7e YE a | 1S ec ia le 
lest layer yield ratio O..978 O. S69 G37) 


Note: the x’s and +’s denote the sSéparate groups of impact 
pointe délineated in Figure 5. 


Theadeation, the angular and energy distributions of the 
150 trajectory simulations were identical in form, and 
nearly sao numerically. Since all the measurable results 
were very Similar, it was deemed prudent to conduct all 
subsequent simulations using only 150 trajectories (those 
Gereagnated by an x in Fig. 9S). 
4. Adatom and Vacancy Placement 

Adatoms are placed at the 3-fold bridge position 
(Fig. 6) of a clean Rhn(1ili) surface by designating, within 
the eimulation, their relative position from a surtace 


layer atom, as measured trom atom center to atom center. 
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This corresponds to 2.1874 in the negative y direction and 
1.5604 in the positive =z direction from the surtace acom. 
For simplicity, an adatom will be designated by the number 
of the surface atom with which it is associated (e.g., 
adatom 108). See Figure 6 for an example of this 
designating system, and Figures 3 and 4 for the atom 
numbering system. 

Vacancies are treated in a much simpler fashion. 


At those positions where a vacancy is desired, the atom is 


simply removed from the lattice. The vacancy location is 
designated by the original atom’s number (e.g., vacancy 
aoe 


A perfect crystal lattice is always assumed in 
these simulations. No attempt is made to account for any 
distortion trom the equilibrium position ot atoms in the 
lattice due to the introduction of adateams or vacanewe—. 
For example, the relaxing ot surrounding atoms inward 


toward a vacancy. 


Ga POTENTIAL FUNCTIONS 
i Backqround 

At present, the knowledge about interatomic 
potential is relatively qood at large separation distances 
Ci eeome, ani the region otf equilibrium separation and 
gQreater), due to the large amount of experimental data 
available. Also, the theory lus Satisfactory at very small 
separations where nuclear scattering is predominant. 
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However, there is a serious lack of information in the 
region between the two extremes, whiten is normally 
compensated for by extrapolating either high energy or 
equibrium potentials based on a reasonable physical model. 
The potentials used in sputtering models lie in this 
intermediate region and are therefore very approximate. 

A typical two-body interatomic potential consists 
of an attractive part at large separation approaching a 
Minimum in the region of equilibrium separation, then 
becoming repulsive and increasing rapidly as the 
interatomic separation decreases further and the closed 
shells of electrons of the atoms begin to overlap. The 
exact form of the interaction in the low energy region will 
differ for two atoms interacting in a crystal and in a 
diatomic molecule since in the former case it is influenced 
by atoms in its neighborhood. 

Like the analytic theories, molecular dynamics 
Simulations treat atom-atom repulsive interaction as two- 
body central forces. The assumption is made that the 
forces and potential energies are functions only of the 
properties of the atoms or ions and the internuclear 
separation. During an MI timestep, several atoms may be 
exerting forces ona single particle; so the resultant is 
calculated as a sum of the individual pair-~wise 


Pier eettone, improvingmehe calculation. 


S57 


Attractive forces are included in this model, even 
though they have little effect on the particle dynamics. 
They do, however, dominate the processes by which atoms 
escape the surface; so their magnitude strongly influence 
the atom yield from a single trajectory. 

The cohesive energy and atom separation of the 
material determine a well depth and "width", establishing a 
reasonable approximation to a pair-poatential well, but an 
attractive force range must be assumed. Prior temehas 


study, test trajectories were run uSing nearest-neighbor, 


NreCal):, interactions, next nearest-neighbor interactions, 
NN(2), eter, source to gin 4 In the fec lattice, dittere jee 
beyond NN(2) are very minor, So normally a potential 
function is truncated between NN(2) and NN(Q). The well- 


depth of the attractive potential function is adjusted so 
that the calculated potential energy of an atom in the 
center of the target equals the heat of atomization of the 
target material. 
Zi Punnett ional Fowm 

Since the early days of™enhe statistical theorirea =a: 
has been predicted that the sputtering yield would be 
relatively insensitive tor ee form ° of “a ~partretera. 


interaction potential —funeeten CRef 


Cs 


: PEW: Sey e The 
critical assumption in any sputtering calculation continues 


to be, however, the choice of the potential functions. 


eRe 


The potential functions used in this simulation are 
empirical in nature. They do, however, provide amore 
Gealbtstic view of the atomic interactions within a cascade 
than potentials derived exclusively from theoretical 
considerations which themselves contain approximations. 
The empirical potentials are based on simple analytic 
expressions which may or may not be Justifiable from 
theory, and which contain one or more parameters adjusted 
to the model or deduced from crystal equilibrium data. 

In general, the atom-atom potentials are compound 
functions consisting of a repulsive "wall", joined smoothly 
to an attractive section, the "well", by a cubic spline 
fime t 1 On. Much of the recent work has been done with the 


Born-Mayer function repulsive potential [{Ref. 57), 
V(r) = A exp (-Br) 


The parameters can be thought of as two degrees of freedom 
in the function corresponding to an ion size and hardness. 
Here the size is the separation at 1.0 eV, and the hardness 
is the pre-exponential factor A. 


Many calculations have also been performed with a 


modified Moliére potential. The original Moliére form is: 
V(r) = ((ZyZpe°/a)/(r/a)) g(r/a) 
men ajeea | O.nso exp(-O. 3 r/a) + 0. 5835 exp(-i.e r/a) 


+ 0.10 exp(-6.0 r/a) } 


a = O. 8853a,/[(Z,1/74 + Zai/ejye/s (Firsov screening length) 
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where a, is the Bohr radius, and Z4 and Zs are he avomuice 
numbers of the ion and atom respectively. Because this 
function is intended as an approximation to a Thomas-Fermi 
potential func Lom, Torrens [Ref. aT] has shown 
justification for using the Pirsov screening seme ta 4, oo 
an adjustable constant. The resulting “2UNClCilOn 22 3cal veer 
moditied Molicre funetion: a’ = Ka. 

The attractive portion of the compound potential 


function Nas the Morse potential tunct liom 109m, 
V(r) =).Dele"P<-25 (2 ae ee a ee eee 


with D, the well depth, ra the equilibrium separaticnyo: sag 
atom pair and a controlling the well width. 
3 Potential Functions Used in the Model 
For this study the atom-atom potential was chosen 
a3 a repulsive, modified Moliére potential joined to an 
attractive Morse potential. The composite pair potentials 
Vig between the ith and jth atoms separated by a distance R 


1s given by: 


Vij = A exp (-BR) Re<uRe 
Vij = Co + CUR] Conse can Ra <R< Ry 
‘a = Delexp<-2@a(R -RS)> = @ exp<-a(R -Re)>) Rp < R < Reo 
vio R SORE 


The functions are Jjoinéa@ smoothly oy eee uotee splines ae 


Polnts kk.) ana Rp: The modified Moliére function is 
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Meencaced at Ro so only NN(1i) SG (cme chactLOons are 
included. 

A Born-Mayer/Morse composite function used by 
Garrison fret. Dole lS ~pmovided —~for “comparison in this 


Stulay, and is of the form: 


Vij = A exp (-BR) Bo 
Wi} = Co “ C,R + CoRe + C3R3 Ra < R < Rp 
ae ‘erelrexp<-ca(R -Ro)> - 2 exp<-a(R -R.,) >] Rp < R < Re 
The interaction of Ar* with Rh is represented in 
both cases Dy? a “wurely reptll sive, unmodified Moliére 
Pume tlOn: 
jus [(Z,Zae°/a)/(r/a)) g(r/a) R < Ra 
Ne =euO R > Rs 
The parameters that define the potential functions 
used in this study are presented in Tables 3 and 4. The 


modified Moliére potential functions for both RH74 and RH7X 


are created by multiplying the Firsov length, a,by 0O. 740. 
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TABEEWs 
POTENTIAL FUNCTION PARAMETERS 


Ries Rhi-kh Rh-Rh Ap RA 
Designation RH74§ RH7X?’ RHBM$ 
A (keV) 0.0 ee, 54.330 One 
Beas to) O7ee) orm) =. O88 0.0 
D. (eV) O. 8237 O27 S29 O26 257 OFO 
R. (A) 2.7350 pany fe) G, Ze Gao 8448, 
Gl, Ange) _ 1.080 1. @86 ie Se) 0.0 
Re uD On 73 Oey S' 0.78 Ze 
ee ee, O. 80 One® Oma 5 Pape G, 
Ra-( US San VAG, 2528 2.40 Zee) 
k 0. 740 O. 740 ae 1a @) 
TABLE 4 
CUBIC SPORE esr Are ter 
Rh-RA solalsde Rhee 
RH74 RH7X RHBM 
Co -4. 264x10? = eo eS 4. 469x103 
Ci i 6eo x110- Ze eee) ~1. S4ex86- 
Ce po ee =2.0 25 10° 1 Je 2xiG; 
C3 9.299% 10> 1235x185 =7.. 1 40xa Oe 


The Rhodium-Rhodium potential functions and forces of 
RH74 and RH7X are compared graphically with the RHBM 
LUNCTION An Frgures 7—- 7e. The force between rhodium atoms 
is determined by taking the derivative of the analytic 
potential functions. The two force functions are then 
joined by a log-linear spline, which causes the distinctive 


break in the graphs as seen in Figures 8 and i10. 


a Function used in 500 eV Ar* simulations. 
cs Function used in the 2 keV Ar* simulations. 


= Garrison’s Born-Mayer/Morse potential function 
[kRet<)-se 


eee Toe 


A. OVERVIEW 

The results presented in this study are the product of 
Poets = Of WOYN simulations totalling nearly 100 in 
number. Two basic avenues ot investigation were pursued to 
determine the effects of smail scale surface defects on 
sputtering. They include: i) the effect of the defect size 
On sputtering; and ii) the effect of the distance between a 
Single surtace defect and the representative area. The 
second area of investigation has provided a measure oft the 
SBurtace defect density required to induce a noticeable 
change in the sputtering. 

The variation in sputtering as a function of defect 
Size was examined by comparing a series of Simulations that 


included a clean Rh{111i} surface and surfaces containing li, 


Zee}, and 7 vacancies or adatoms. The locations of the 
defects, relative to the impact '§ area, are depicted in 
Figures ii and iz. Deviation in sputtering was similarly 


studied by comparing the results from a clean Rh{1iil} 
Surface with simulations that included a single adatom or 
vacancy at various distances from the impact area. Both 
sequences of simulations were carried out at bombarding Ar’ 


energies of 500 eV and 2 keV. The investigation of several 
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Side issues such as the comparison of potential functions 
and containment criteria (see Fig. 1) were also conducted. 


Voluminous amounts fad data were collected trom this 


series of computer simulations. Presented here are 
portions of thie antormatton, assembled into a series of 
tractable comparisons, inc bucaag : ee, the effect of 


# 


vacancies as compared to adatoms, and 11) the ditterences 
between simulations run with ion energies of 500 eéV and 2 
keV. The comparisons are made amonget various 
experimentally measurable quantities such as the energy and 
angular Gi strepur Len ot sputtered atoms to aid in 
discerning trends in the data. Various mechanisms are also 
discussed to account for the variations observed in these 
quantities. 

Table 5S provides a summary of the operating parameters 
used in the majority oft the Rh{111i}/Ar* simulations 
conducted tor this study. 


TASBE SS 
QDYN OPERATING PARAMETERS 


Ion energy (eV) 500 2000 
Ion angle of incidence 2O¢ = | 8 
Lattice dimensions Leases 12 20 XO sa1e7 
Rh-Kh potential vtunce ven RH74 RH7XK 
Ar-RAS potential fune@ren*? Born =-Mayer Born-~Mayer 
Number of trajectories Ga) 8, 150 
Representative area 

pin point Catom +> ea 108 


Modified Moliere/Morse composite potential function. 
See Tables 3 and 4 for function parameters. 


1° See Table 3. 


4-4 


Ee SrultTerinGg YIELDS 
Ties oT or Containment af the sauttering process 


w2thin the microcrystal otf a simulation, Harrison has shown 


that: “Only the First few layers af the target are needed 
tae conpletely explain ALL experimental souttering rasults, 
Pep os 7 ee Wee tL cD re or as Mae ERet: So) Since absolute 


containment 15 not achieved in this model the total yield 
welt be Low when compared to experimental results. All 
other computed quantities, however, will be very close to 
Pafeeeow tah inite Crystal" values for the lattice sizes used 
in these simulations. Relative yields are therefore 
reported, with the values being normalized to the clean 
surface results. 

Figures 13 - i6 detail the effect of a single defect on 
the normalized yield as measured by the defect’s distance 
trom the center of the representative area. The center of 
the representative area was chosen as the reference point 
Cawieme than the pin point, since the impact points are 
almost uniformly distributed about the center. 

For the low energy case a signiticant portion of the 
yield is initiated by thé recoiling primary knock-on atom 
Gara Gr "target atom". The comparison at Ar* energy of 
200 eV trig. 13) indicates that a vacancy has a greater 
effect on the yield than an adatom does. This is to be 
expected, eince the incoming ion has the opportunity to 


transfer its energy to a second layer atom through direct 
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cod ilistenm, thereby involving atoms much deeper in the 
lattice. Subsequent collision sequences would distribute 
the energy through layers deeper in the microcrystal than 
would normally be expected. The likelihood of a surface 
atom having enough energy imparted to it to eject it from 
the surtace (about 25 eV) is thereby decreased. 

The presence ot a peak at 1.5 LU tor the adatom curve 
at Ar* enerqy of 300 eV (Figs. is and 16), correspencing =o 
an adatom at positions 42 and 61'!', cannot be reconciled by 


a Similar mechanism and further study of the detailed 


tf) 


collision sequence is required. 

At Ar* energy of 2 Bey 3 owe Signiticant yield wre 
attained from both the recoiling target atom and the 
scattered ion, where the "ion" yield is increaced by 


contributions from second and third PKA’s. Again, the 


E 


(3) 


moa 


i) 


vacancy has a greater effect on yield. This (ef 2ec t. ma 
probably caused by the same mechanism described above, but 
in this case the energy is lost trom the microcrystal by 
having more atoms ejected from the sides and bottom of the 
microcrystal, as compared to a clean surface. 

The ability of’a single adatom Cat, pos it tenes LO Se 
increase the yield can be explained by visualizing the 
adatom being driven into the surface at an angle by the 


impactaag loam The propagating secondary collisions would 


See Fig. 3 for adatom positiona. 
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Pema to be closer to the surface, than an ion impacting a 
clean surface. The recoiling atoms would then have an 
opportunity to transfer a greater portion of their momenta 
Bemeene Surface atoms directly, causing more atoms to be 
ejected. This effect is not present at S500 eV since the 


Pero More Likely to be reflected off the adatom. A&A briet 


ee 


comparison ot reflected Ar* 1S presented ain Table to 


Support this hypothesis. 


TAB GE 6 
NeiPoek OF REP LECTED Ars ONS 


Ton energy (eV) 300 2000 
Clean surface 53 46 
1 adatom* 96 36 


*Note: adatoms are at the same positions relative to 

the representative area. 

The scattering of data points for the 2 keV adatom 8 
Crees. Lt and i6) Was investigated in some detail. The 
data points at approximately 1.8 LU correspond to adatoms 
Mmegociltions 9D (normalized yield of 1.028) and i20. While 


these adatoms are the same distance from the center of the 


representative area, they are not the same distance from 
each of the individual impact points. Where there existed 


a Significant difference in yield per single ion, the 
sequence of collision was’ traced. The most Significant 
ditference between the two trajectories was the position on 
the eurtface where the ion actually made contact. This was 


caused by a slight difference in distance between the ion’s 
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starting position and the adatom. The lang ragge oft? che 
Ar-Rh potential function (2.2 LU for cutoft) .eauseag. 7 
incoming ions to follow difterent paths toa impact whien 
yielded two different collision sequences. This was 
confirmed by reducing the Ar-Khee@arential! Sflunet vongerueanw. 
to 1.6 LU, at which time the tragjectery resv bee teeter 
exactly. ‘No trend can be observed for the dispersed data 
pointe as a whole concerning their relative positions toa 
one another. 

In general the qraphs asymptotically approach a 


fect 


IT 


normalized@yzeta” et ita Thi2] cGenressond=Stesethe ad 
being moved tarther from the area of collision events until 
it has no effect on the trajectories at all. The distance 
at which the curve reaches +5 percent of the clean surface 
results is an indication of the density of defects required 
to visibly attect the yield. This takes into account an 
average error of 324 as measured by the deviation of the 


curves from 1.0, and uses only the upper data points of the 


<“ keV adatom curve. Table 7 presents the approximate 
distances as determined from the Stratqae line 


extrapolation of Figures 13 and 14, and the surface density 
corresponding to a Gircle %emeleseine a Geteet "at ensae 


pesitioan. 
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TABLE 7 
SURFACE DENSITY OF DEFECTS 
CAUSING +5% DEVIATION IN YIELD 


Distance of defect Suriace Penei cy 
from center of (defects/iA ?) 
representative area (LU) 


Se@neV, 1 vacancy Loo Oe ta3 
200 ev, 1 adatom Ooa7 aes a5 
2 keV, 1 vacancy dee EO a7 
2 keV, 1 adatom OAS se 0.0492 


To add some credence to the predicted defect density 
pestlts, Simulations were run at 2 keV that evenly 
distributed the defects over the surface. Figures 17 and 
18 depict the positions of the vacancies and adatoms 
respectively. The closest defect to the center of the 
representative area is actually outside the distance listed 
My Table 7 in order to provide a check an the outer limit 
of results. For the vacancies a raw yield of 957 was 
obtained, versus 1022 for a clean surface; a 6.4% decrease 
mothe yield. The adatoms produced a raw yield of 1001, 
wen 15 only a 2.1% decrease in yield; not quite the S% 
that was desired. 

The difficulty with the adatom case is that the next 
closest position, nearer the representative area, to place 
an adatom would have located it within the representative 
area. This was felt to be too biased a test to confirm the 
cesults, since it would have essentially duplicated tne 
Perens t Simulation, and was therefore not conducted. 


Dezpite the lack of confirmation of the adatom results, I 
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believe the the defect densities necessary to produce a 5% 
change in the yield, as Shown ia Table 7, co be Seneaceureace 


estimate. 


an ATOMS EJECTED PER SINGER er 
With the advent of computer modelling the ability to 


determine a distribution of atoms ejected per single ison 


1a 


(ASI) became available. This quantity cannot be determined 
experimentally but it does present some insight into the 
mechanics of the system. 

Fiqures 132 and 20 present a series of comparisons 
Within a Single ion energy regime tor: a clean surtace; 1 
and 3 vacancies; and 1 and 3 adatoms. For the S00 eV case, 
little difference is seen between the clean surface 
distribution and both of the adatom results. This is as 
expected trom the sputtering yield results, where adatoms 
have been shown to have little impact on the total yield. 

Again, as expected from the sputtering yield results, 
vacancies have a Siqgnitxrveane etfect on the ASI 
distribution. The mass of the distribution tends to shift 
to the left as the size of the pit formed by the vacancies 
gets larger, With the low yield trajectories becoming 
dominant. 

At 2 keV much the same behavior can be seen happening 
to the. Vadtstributwvens, With the same effect now being 
observable for the adatoms, but to a smaller extent when 
compared with the vacancy cases. 


Se. 


am alternate form of displayimg the ASI distribution 
Which accentuates the atom ejection process is shown in 
hegure 21. eS ainiser of rings ~im each cirele is 
proportional to the ASI value for an impact point at the 
center of the circle. The resulting pattern indicates the 
areas of the representative area which produce high yield, 
eueeeom, by the relative darkness of the circle. Since the 
RH{1ii} surface does not contain any true SIDS Vygis esters (= 
is no significant difference noted in any atom wield per 
Petes Eaint diagrams created during this study, Figure 21 


being provided as an example only. 


BP Peaer lon TIMES 

Pig@iees 22 and 23 provide a comparison of ejection time 
Sicoerisucions for Rh, within an ion energy regime, as they 
are affected by 1 and 3 vacancies, and 1 and 3 adatoms. At 
[900 eV there is a slight shift to lower ejection times for 
an increasing number ot vacancies, and a very considerable 
shift for the adatom cases. The shift to lower ejection 
times at 2 keV is barely discernable in all cases. 

The most significant fact that can be seen from these 
Mee cicns is that more than 95% of the atoms ejected 
during the simulation occur prior to 250 femptoseconds 


@e ts = 10°-'> sec) for the 3500 eV cases and prior to 300 fs 


at 2 keV. 


Sil 


Ew LAYER ,CUEeP DS paras 

The layer yield ratio is defined as L, = N./N, where Ma 
is the number of ejected atoms which originally resided in 
the mth layer of the target, and N is the total number ot 
atoms ejected. Lex values at SOO eV and 2 keV are plotted 
in. Pigurses” 24 Sane = At both energies, the adatoms 


"new" first layer causing the yield 


( 
fb 


eftectively act a 
ratio to decrease more rapicty shan the vacancy eestor. 

For all cases, except at low energy and high number of 
defects, the ejected atoms come primarily from the tirst 
target Laver. This agrees with Harrisons statement that: 
"Under all conditions 85% or more of the ejected atoms come 
trom the surtace layer of the target, and essentailly all 
of the remainder come from the second layer". ([CRef. 595] 
However, as the number of defects increase in the low 
energy case, the target is more able to accomplish a 


momenta reversal and have atoms ejected from deeper layers. 


ae RJECTED J ATON ENERGY sbi st rious 

The ejected atom energy distributions were created by 
tallying the energy of all Rh atoms sputtered from the 
surface of the target (whether the Rh was single astems or 
multimers) into energy bins O.5 eV wide. Pigugese26)--— 24 
are plots of the indivaeual eae An attempt to plot a 
smooth curve through the Bein points is not made, since 


Significant noise is present in the plotted values. 
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Figures 26 and 27 present information obtained from 2 
keV simulations. It can be seen that there is little 


difference between the energy distributions as the number 


of defects increase. These distributions, as a group, can 
be viewed as one ensemble. By averaging the energy values 
moa each bin, their statistical plausibility could be 


improved. 

Pigure 28 as compiled fram simulations run at S00 eV 
With GOO trajectories each, since those completed using 150 
trajectories did not contain enough information to discern 
Snything beyond the moise. 

Everotlowing the data points to be interpelated by eve, 
these energy distributions can be seen toe Gfollaw, > 2 
Maxwellian Lorm, with a Maxima at about 2.5 eV. 
Additionally, there is a slight secondary peak at 9.3 eV 
tee the 2 keY cases and one at 8.0 eV for the 300 eV graph. 


At this time no explanation is ftorwarded to account for 


these secondary peaks. 


Eh, EeeereD ATOM ANGULAR DISTRIBUTIONS 

Semee 4 major thrust of Garrison’s work (Ref. 38] has 
been to explain the location of the peak otf the polar angle 
distributions, aS compared to experimental results (Refs. 
48 - 50], an extensive group of plots (Figures 29 - 41) has 
been developed to compare the etfect of i0On energy and 


Seteat oitze on the distributions. 


oS 


The polar angle distributions are created by tallying 
only those ejected Rh atoms with energies greater than 5.0 
eV and summing over all azimuthal angles. Diserinana = 1G 
against the very low energy species will help align the 
results more closely with the experimental measuremencs. 

At 2 keV the shiitt in the @ietributitoen is verye= lo 
With the vacancies forcing the distribution to smaller 
angles and the adatoms to larger ones. This effect ads nuct 
more pronounced in the SOO eV results, with Figures 37, 40 
and 41 showing the shifts most clearly. In general, as the 
number of detects increase; tram a clean Suriece theetan 


2, 3, and 7 defects'*; the amow@nt of the shitt in the peak 


He SPOT PATZIERNS 

An alternate method tor displaying angular distribution 
is the sputtered atom spot pattern. Unlike the angular 
Gistribuwetoans discussed en section G above, this method 
reveals both the azimuthal and polar angle dependence of 
the sputtered atoms simultaneously. 

Little difference was noted between the spot patterns 
of the simulations conducted for this study, theretore only 


two examples are presented here for illustration (Figure 


we See Figs. 141 and” 12 tor) “erative esesitloneen 
defects. 


24 


2) Of note is the expected six-fold symmetry of the 
Pomeern Characteristic of a fee(ili). 

These graphs were produced by plotting the ejection 
angles otf all sputtered atoms. To Obtaanispet patterns 
more Similar to those experimentally obtained, such as SIMS 
data, Many more trajectories would be required (on the 
Seager of 1000). Atoms ejected with energy less than 5.0 eV 


would then be eliminated toa produce a sharper image. 


ig Pee enArem FOrENTIAL FUNCTION 


fas Majority oT theoretical WeoCrmgen tne Pri iit paoas 


Ui 
ae 
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has been conducted at Pennsylvania State University 

By 2B. Je Garrisenm ([Rexr. See ti .congumrian wich the 

experimental studies of N. WinoGrac, et. al. C[Refs. 48- 

=) On - Due to the close collaboration between Professor 

meee teson and Professor Garrison, it was considered de 

riguéur td compare the results obtained using the Born- 

MWeyes Ra-Rh potential function ‘(RHBM) with the modified 

terre re potential function, which is used throughout this 
Study (refer toa Chapter 3, section D.3 for details of thes=s 

ue r loans). 


Simulations were run on clean Surfaces only, at 10n 


energies of S00 eV and 2 keV. Significant differences in 
results are apparent. Summaries of the graphical 
comparisons are listed below, with the statements being 


couched as to the effect of the RHBM function when compared 


to the Moliére function. 


Cn 
Cn 


i) 


2) 


3)) 


4) 


Figure 43. Atoms sputtered per single ion (ASI): 
At 500 eV a noticeable shift to lower ASI is 
observed, with only a @miner shaiit apreartis toe 
kev. 

Figure 44. Atom ejection time” “distr ilaweron. At 
both energies a shift ta lower ejection inet 
evident, with the times “required tor "so, Sjeeee 
being reduced by 60 fs. 


Figure 45. Ejected atom energy distribution: The 


data presented has a low order of confidence due to 


insutticient atoms ejected and the noise presen =e 
the “date errpurtiens: However, some details can be 
determined. There is a slight shift to lower 


energies for both S00 eV and 2 keV caseg, with a 


corresponding reduction in the maxima. 


Figures 46 £4and 47. Ejected atom angular 
Si St Yao c Lemon. These tigures present the most 
dramatic difterences between the potential 


functions, with the RHBM results being very sharply 
peaked at 25° for both energy regimes, compared to 
broad distributions peaking at larger angles for 
the Moliere results. 

Figures 48 and 49. Atom yield per impact podnee 
These figures are presented IGE a general 


information. No significant differences were noted. 


Table 8 presents the raw yields for each of the four 
Simulations compared, and a percent change in the yields, 
to emphasize the magnitude of the effect the choice of 


Peeential function can have. 


TABLE 36 
COueenr lOO OF SkAW Yiteps 


Ion Energy (eV) 1G) 2000 
Yceid (RHBM) e765 POs82 
Yield (RH74, RH7X) 457 1022 
Percent change =18 ve IL 


a4 


Whe CONC USUaNS 


While this study has madé no attempt to compare the 


simulation results to experimental results that are 
currently available, it has demonstrated that the presence 
of defects  Gmy single eh ghz guns Ul Surtaces can preduce 


measurable difterences in the quantities that characterize 
Sputtering. By investigating yet another variable in the 
Simulation of sputtering, this study has laid to rest some 
of the criticisms leveled against the simulations in 
general, e.g., a Ssimulation’s lack of accuracy due to the 
incomplete modelling of the system. In the future, the 
placement of vacancies and adatoms may be considered when 
any system is modelled in order to provide: i) a more 
accurate prediction of experimental results, ii) possibly 
greater correlation with the existing @xperimental results, 
and “i asam) greater insight into the actual sputtering 
mechanisms as simulations more closely reflect actual 
sputtering events. These statements are predicated upon 
the assumption that future experiments will be conducted 
with the precision evident in the EARN experiments 
conducted By N. Winograd, “st. al. [kets 44)-s 50. 

The correct choice of a moditied Moliére/Morse potential 
over a Born-Mayer/Morse potential was verified when the 


results of the simulations were compared to preliminary 
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BARNe results (Ref. S&l. Garrison’s simulations, using the 
RHBM potential AU eo tot) COmeE@e rect ayn, edict the 
peaks occuring in the angular and energy distributions, 
while the simulation results of this study are a very good 
Seeasoltatrvanm tor the angular distributien only. Correlation 
of the energy distribution is left for later studies. 

Pree noTe sot soetemtc tal functien does, however, cloud 
the issue of the etfect of the surface defects. The choice 
of an appropriate potential function is still considered to 
beware Glack art", Sinmee little is known of the behavior of 
potential functions in géneral at the low energies at which 
the sputtering processes operate. The choice of an 
appropriate Aeon only has to fit a few parameters, such 
as the surface atomization energy of the crystal. Beyond 
that, the person creating the model’s environment has wide 
latitude. 

The magnitude of the effect of the potential function 


choicé is at least as great as that of the surtace detects 


reported here, so it can be argquéd that possibly ar 
incorrect potential function was chosén tor this model. if 


da not believe that this is born out by the results 
Seecsentead in this study. Surface detects exist on all 
surfaces used in experiments and simulations should take 
[ime ato acceunt. 

Future studies into the Rh{1lli}/Ar* system using QDYN 


Simulations should increase the ion energy to 5 keV to 


So 


match the experiments currently being conducted in this 
area. Also, more attention should be paid to investigating 
the mechanisms that may be responsible for the results 


observed, by tracing the progress of the collision events. 
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APPENDIX: 


ELGURES 


This appendix includes all figures referred to in the 


peayeet Chis report. The 


are the measurement of 


rhodium atoms. The energies 


refer to the impacting 


Simulation. 


data presented in these 


Various parameters of the 


argon ion 
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Fig. 6 Adatom Position. 


val. 





Ve 
3 


POTENTIAL ENERG 


RHODIUM—RHODIUM POTENTIALS 


a olialA 


0.20 0.90 0.75 1.00 120 Beat ise 


7 





SEPARATION (LU) 


Rhodium-Rnodium Potential Function. 


U2 


0.000001 


0.000000001 


RHODIUM—RHODIUM FORCES 





Oo 0.50 0.75 100 L.2o 
SEPARATION (LU) 


Fig. 8 Rhodium-Rhodium Forcea. 


73 


RHODIUM—-RHODIUM POTENTIALS 


~ 


Da 


POTENTIAL ENERGY ( 





0.20 0.50 6.735 1.00025 ete ome ooo! 
SEPARATION (LU) 


Fig. 9 Rhodium -Rhecitum Potted (a tee eed. 


74 


O.000001 


0.0000001 


‘P ‘ 
— -_— 
= 
—) 
7 
a 
(eae 
| on, 
one 
SS ve 
f —- Caos 
——| e 
NS 
Nene 
al 
= 
= 
al 
oe, 
ee 
an, 
— 
—_ 
wee! 
a) 
re 
ae 
— 
—) 
e 
a) 


RHODIUM—-RHODIUM FORCES 





Owes e, 0.50 Oi 1.00 ies 
SEPARATION (LU) 


Fig ss lLOeeRnOGtUune=Rnocium horces. 


TS 


t.50 


i VaCaie, 2 Vacancies 


+ 


3 Vacancies 7 VaGanerves 


Pig. 11 Relative Vacancy Positiona. 





Vs 


; 1 Adatom _2 Adatoms 





| | 





Cf 


_/ Adatoms 





Fig. 12 Relative Adatom Recae dens . 


TF 


OS 


(V1) JONVLSIO 
O'v ove OZ 


Ad 009 ‘ADYINS HV “WOlvav | 9. 
Ad 00S -ADYANS HV ADNVOVAL oO 
QN3941 


VAYV SALVINASSYddY AO YSLNAD WH 


ADNVLSIG SA CTsIA 


Git 





80 


6'0 


om 
Q1alA GAZMVWYON 


Zt 


Yield versus Distance from 


1S 


Pasa 


Center of Representative Area. 


78 


In 


(M1) JONVLSIG 
O'r Oe Oz 


{ 


AA C -ADYANA HV “WOIVOUV I oO 
Ady ¢ -ADYUSNA HV CAONVOVAL © 
QN3H941 


VIUV JALIWINASSHd3Y 4O WIND WH 
AINSI “S/\ GTA 


Ol 


0'O 





60 80 L'0 


Q1slA G3SZIIVWYON 


a 


ot 


Yield versue Dietance from 


Center of Representative Area. 


14 


Pac, 


a=, 


O'S 





(N1) JONVLSIO 
O'P Ove oz Ol 


U 





Ad 00S -ADYANS YV CADNVOVAL oO 
INES 1 


VINV JALLWLNISSYdaY JO YIINIO WS 
ADNVLSIGC SA GAIA 


O'0 


O'1 
QialA G3aZllvWyON 


ene, 


8°O 


60 


Lt 


Zl 


Yield vereuse Digetance from 


es 
Center of Representative Area. 


hoon 


80 


O'S 





(11) JONVLSIG 
O'r O€ Oz 


AIX Z:ADYANT HV ‘WOLVGV Lo 
AJ 00S :ADYUSNA HV ‘WOLVAV Lo 
GNa971 


VIHV JALLWINISTHdaY 4O YAINTO Wd 
FO] NAvAbs] (@ mesy AMG MEI PAN 


O'l 


00 


8°Q LO 


60 


O'l 
Q1glA GQaZlIVWYON 


Lt 


ct 


Yield vereus Distance trom 


16 
Center of Representative Area. 


Pala: 


81 


——— gy ——~ a > = = a ae 


emcees ey, ——wy —— | 
—— \ — A ——— aaah 5 = 
oN oe Oa ee — 
; Ve Y : , y F f if ; oe 
\ A A A { A { | ‘ i { aN 





(ik NX i XD oe 





Pig. 17 Distyatbwtec ly acanie ae aee no. 


82 





Fig. 18 Distributed Adatom Positions. 
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Fig. 19 Atoms per Single Ian Distribution. 


84 


ATOMS PER SINGLE ION, AR ENERGY: <KEV 


XQ 

- | 2 ADATOMS 
= A Wor. 17m 

. 

5 Ly 1 ADATOM 





3 VACANCIES 


FRACTION 


GEO (Ue aU ye UU Oa Wc eee Ol Ole 


1 VACANCY 





GneaN SURPACE 
LAE til WITT. A 


3 10 ie) 20 29 
ATOMS Pee S UN CIEE, Jd omy 


Fig. 20 Atome per Single Ion Distribution. 


85 










| } () oe oe oe OO ye oa on 
a7 ©) Gop 1909 (©) ()) \O OOD E (Oo) (ce (on) &) VO) 4O)) 
) Hoy) G (MO) (on) ( >) ©) 1Or (OME >») Or ‘y)) (« neo) 


>A 


ve ( ey y (ous ( { (i) (ls iy) ( (i: a) | Gs ai) Kony | Kony (ec i») ( fry) ( )| (on veony ( oy C Fem) (@a} Goo Won 


ee ‘é) TeX (© 


° 


(091) KG) CHO) € Y 


— 


— 


‘ey ©: a) V €. )) 1 fe.) ) fy) O) & 1o) ( (Cori) (( | (& D a (o) (ws lj G )) Cy Kon ot D 


igs Ss Z- MN - tH 


(OO) NO} (« nh (O) oO) ) (ONG ©) (onc O) C \( (9) ©) ( op Oe WH Cy ( () ey) 


~ 2) SSE 


(6; PNY OMG (« %) O} (o) (0) CO) ie ) ee) roy (om © (2) (©) (O)) (( (oo ( i 3) 


~ ok 


IO VOQVOOOOOO o 6) oye DOO 


T<- Bo 


es ()) (c heyy ( (° MC O) (Ne »») ©) | (Colic) )) (©) (oo) Cre (@) () (@) (eo) 


a 


Qo ) Geom (OV ( 3) Oe oy) (ox (wn) OM ©) ory ( an On ©) COG (a a Ot 







Atom Yield per imple emt odmus 








eee 4 Ne r©& i) ye aN ()) (9) (( () Grex ) fe) or 
5) ( 
Oo 


O)) )\e Dy) (co) Ke ©@)(O) (@) QO OO € 20) ©) 
Cr1©) i reac ee ant ( \(c ©) ®) 
VOCE 6) OOOO ) Oo 





ATOM YIELD PER AWiPACT PORN 
et 


Pig: 


500EV RH(111)/AR<111> CLEAN SURFACE 


86 


TIME DISTRIBUTION, AR ENERGY: SOOEV 


HRACTION 


OMS O “Oc (0°02 Oe em e ee Onl 2 


3 ADATOMS 





Ne 


1 ADATOM 


1 2e|| a S 


3 VACANCIES 





1 VACANCY 


CEEaN SURFACE 





200 300 400 
EJECTION TIME (FS) 


Fig. 22 Atom Ejection Time Distribution. 
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